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The 6 Li abundance observed in metal poor halo stars exhibits a plateau similar to that for 7 Li 
suggesting a primordial origin. However, the observed abundance of 6 Li is a factor of 10 3 larger and 
that of 7 Li is a factor of 3 lower than the abundances predicted in the standard big bang when the 
baryon-to-photon ratio is fixed by WMAP. Here we show that both of these abundance anomalies can 
be explained by the existence of a long-lived massive, negatively-charged leptonic particle during 
nucleosynthesis. Such particles would capture onto the synthesized nuclei thereby reducing the 
reaction Coulomb barriers and opening new transfer reaction possibilities, and catalyzing a second 
round of big bang nucleosynthesis. This novel solution to both of the Li problems can be achieved 
with or without the additional effects of stellar destruction. 

PACS numbers: 26.35.+C, 95.35.+d, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 

It has recently been pointed out (e.g., [l|) that the 
abundances of both 6 Li and 7 Li observed in metal poor 
halo stars (MPHSs) are not in agreement with those pre- 
dicted from standard big-bang nucleosynthesis (BBN). 
Specifically, the 6 Li abundance as a function of metal- 
licity exhibits a plateau similar to that for 7 Li in very 
metal-poor stars, suggesting a primordial origin for both 
isotopes. This 6 Li abundance plateau, however, is a fac- 
tor of ~ 10 3 larger than that predicted by BBN. A less 
severe problem exists for 7 Li; the BBN value based upon 
the baryon-to-photon ratio fixed by the WMAP analy- 
sis 0] of the cosmic microwave background is roughly a 
factor of three higher than is observed. 

Moreover, a longstanding effort in cosmology has been 
the search for evidence of unstable particles that might 
have existed in the early universe. It is thus natural to 
ask whether the two lithium abundance anomalies might 
be a manifestation of the existence of such a particle. 
In this context a number of possible solutions to the 6 Li 
problem have been proposed which relate the Li anoma- 
lies to the possible existence of unstable particles in the 
early universe J3|. This has been extended in several 
recent studies [J H, H, 0, H, [§] to consider heavy nega- 
tively charged decaying particles that modify BBN, but 
in rather different ways. In these latter studies, the heavy 
particles, here denoted as AT - , bind to the nuclei pro- 
duced in BBN to form X-nuclei. The massive X~ par- 
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tides would be bound in orbits with radii comparable to 
those of normal nuclei. Hence, they would reduce the 
reaction Coulomb barriers thereby enhancing the ther- 
monuclear reaction rates and extending the duration of 
BBN to lower temperatures. 

Pospelov [i| suggested that a large enhancement of 
the 6 Li abundance could result from a transfer reac- 
tion involving an X" bound to 4 He (denoted 4 Hex), 
i.e. 4 Hejs:(<i,A~) 6 Li. Although this was an intriguing 
idea, Hamaguchi et al. [9( have recently pointed out via 
a more complete quantum mechanical calculation that 
Pospelov's estimate for the 4 Hex(G?,A~) 6 Li cross section 
was too large; this leads to too high a 6 Li abundance. 
Cyburt et al. [1| further motivated this hypothesis by 
identifying the A - as the supersymmetric partner of the 
tau lepton, i.e., a stau, and considered A~ transfer re- 
actions for 7 Li and 7 Be production, too. Although their 
calculation is based on a fully dynamical treatment of 
the recombination of nuclei with X~ particles and also 
of BBN, they used cross sections for all X~ transfer reac- 
tions involving 6 Li, 7 Li, and 7 Be that were too large, as 
we discuss below. Therefore, the calculated abundances 
are to be viewed with caution. Kaplinghat and Rajara- 
man Q observed that the decay of an X~ when bound 
to 4 He would occasionally knock out a proton or neutron 
to produce 3 He or 3 H, thereby enhancing their abun- 
dances and the abundance of 6 Li through reactions with 
other primordial 4 He nuclei at higher energies. Kohri and 
Takayama 0] studied the recombination of nuclei with 
A~ particles, and suggested the possibility of solving the 
7 Li problem. However, they did not carry out dynamical 
calculations involving recombination processes and BBN 
simultaneously. This forced them to introduce artificial 
parameters for the fractions of the captured nuclei, which 
turn out to be different from the fractions obtained by 
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solving the recombination plus BBN fully dynamically. A 
new resonant reaction 7 Bex(p,7) 8 Bx has recently been 
proposed by Bird et al. [|[ that destroys 7 Bex through an 
atomic excited state of 8 Bx, and the present study identi- 
fies another effect in this reaction that might also destroy 
7 Be. Thus there has been a great deal of recent progress 
in X~ catalyzed BBN in three important aspects: the si- 
multaneous description of recombination and ionization 
processes of X~ particles with nuclei in the description 
of BBN, use of updated reaction rates involving the X- 
nuclei, and inclusion of new resonant processes by which 
7 Be is destroyed. No previous calculation has involved 
all of these effects in a single dynamical calculation of 
BBN in order to study their effects on the 6 Li and 7 Li 
problems. 

In this Article, we present the results of a thorough 
dynamic analysis of the effects of X~ particles on BBN. 
The important difference from previous works is, firstly, 
that we carried out a fully dynamical BBN calculation by 
taking account of the recombination and ionization pro- 
cesses of A - particles by normal and A-nuclei as well as 
the many possible nuclear reactions among them. Sec- 
ondly, the reaction rates on normal and A-nuclei used 
in the present study are based on quantum mechanical 
calculations of the cross section like those of Hamaguchi 
et al. Q, which we believe to be correct. Thirdly, we 
have not only included the important Be destruction 
mechanism identified by Bird et al. [8(, but have iden- 
tified another potentially important destruction mecha- 
nism involving the reaction channel 7 Bcx+p — > 8 B*(1 + , 
0.770 MeV) x -^ 8 B x +7 which has the potential to de- 
stroy 7 Bex via capture through the f + nuclear excited 
state of 8 B. We show that when all these effects are in- 
cluded, the single hypothesis of the existence of the X~ 
particle provides a remarkable solution to both the 6 Li 
and 7 Li abundance anomalies. We can then use this con- 
straint to place interesting limits on the A~ relic abun- 
dance and its decay lifetime and mass. 



II. MODEL 

We assume that the A~ particle is leptonic and of spin 
0, consistent with its identification as the supersymmet- 
ric partner of a lepton. The A~ would be thermally pro- 
duced at an earlier epoch together with A + . Their small 
annihilation cross section allows a significant abundance 
to survive to the BBN epoch. The mass and decay life- 
time of the A~ is ultimately constrained by WMAP and 
the present BBN study. Only the X~ can bind to nuclei 
and the A + remains inert during BBN. The binding en- 
ergies and the eigenstate wave functions of the A-nuclei 
were calculated by assuming uniform finite-size charge 
distributions of radii r$ — 1.2A 1 / 3 fm for nuclear mass 
number A jfO|. When the A~ abundance is very high, 
some nuclei can bind two X~ particles, such as 3 Hcxx 
and 4 Hexx- In that case their binding energies were cal- 
culated using a variational calculation with a trial wave 



function for A-nuclides bound to one A particle, anal- 
ogous to the case of the H^~ ion. 

Thermonuclear reaction rates (TRRs) for all reactions 
that might take place in X~ catalyzed BBN, including 
the A~ transfer reaction suggested in [1] and A~ decay, 
were added to the BBN network code. See Kusakabe et 
al. |ll| for details on the calculations. These were cor- 



rected for the modified nuclear charges and the effective 
mass resulting from the binding of one or two A~ par- 
ticle (s). If the X~ decayed at some later stage, they 
would be expected to destroy some fraction of the nu- 
clei to which they had become bound during BBN. How- 
ever, that fraction would be small We found that 
the inclusion of the A-nuclei 8 Bex and s Bexx (both are 
bound) results in a leakage of the nuclear reaction flow 
out of the light nuclei or A-nuclei to produce slightly 
heavier A > 9 nuclei. This might be an additional BBN 
signature resulting from binding the A~ particles. We 
determined most thermonuclear reaction rates involving 
the A-nuclei by taking account of the lowered Coulomb 
barriers and modified reduced masses. However, as dis- 
cussed below, there arc a number of reactions that require 
careful additional considerations. 

As noted by Pospelov reactions in which an I" 
particle is transferred can be very important in circum- 
venting some normally inhibited reactions, especially the 
4 Hex(c?,A~) 6 Li reaction. Its rate could be orders of mag- 
nitude larger than that of the 4 He(e?,7) 6 Li reaction, which 
is suppressed due to its occurrence through an electric 
quadrupole transition. Hamaguchi et al. have recently 
carried out a theoretical calculation of the cross sec- 
tion for 4 Hex(d,A~) 6 Li in a quantum three-body model. 
Their value was about an order of magnitude smaller 
than that of [|| . This difference can be attributed to the 
use of an exact treatment of quantum tunneling and a 
better nuclear potential. We, therefore, adopt the result 
of [9j in the present study. 

Cyburt et al. {B} estimated astrophysical S-factors for 
the 4 He x (i,A-) 7 Li, 4 He x ( 3 He,A-) 7 Be, 6 U x (p,X~) 7 Be 
and other reactions by applying a scaling relation 
Sx/Sj cx pfao/(w 7 ao) 2A+1 . Here, Sx and 5 7 are the 
S-factors for the A~ transfer and radiative processes, re- 
spectively, cio is the A~ Bohr radius of 4 Hex or 6 Lix, 
Pi is the linear momentum of the outgoing 7 Li or 7 Be in 
the A - transfer reactions, and w 7 is the energy of the 
emitted A = 1 (electric dipole) photon in the radiative 
capture. However, the reaction dynamics are important 
to these results. 

4 He, 6,7 Li, and 7 Be occupy an s-wave orbit around 
the X~ particle (assuming the X~ particle to be much 
heavier than these nuclei). The 6 Li nucleus is an a+d 
cluster system in a relative s-wave orbit, while the 
A = 7 nuclei are a+t and a+ 3 He cluster systems in 
relative p-wave orbits. This difference in the orbital 
angular momentum will produce a critical difference 
in the reaction dynamics between the 4 Hex(d,A~) 6 Li 
and the 4 He x (i,A~) 7 Li, 4 He x ( 3 He,A-) 7 Be, and 
6 Lr>f (_p,A~) 7 Be reactions. Specifically, the latter three 
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reactions must involve A/=l angular momentum trans- 
fer. In order to conserve total angular momentum the 
outgoing 7 Li and 7 Be in the final state must therefore 
occupy a scattering p-wave orbit from the X~ particle, 
leading to a large hindrance of the overlap matrix for the 
X~ transfer processes. Thus, a realistic quantum me- 
chanical calculation results in much smaller S^-factors 
than those estimated in [f|. Therefore, in the present 
study, the above three reaction processes were found to 
be negligible and were therefore omitted. 

Bird et al. suggested that the 7 Bex(p,7) 8 Bx res- 
onant reaction could destroy 7 Bex through an atomic 
excited state of 8 Bx ■ They also proposed that a charged 
weak-boson exchange reaction 7 Bex — > 7 Li+X° followed 
by 7 Li (p, a) 4 He could destroy A — 7 nuclides. We in- 
cluded only the former resonant reaction in the present 
study, although we confirmed their assertion on the weak 
process as will be discussed in a separate paper. 

In our exhaustive study of additional processes re- 
lated to 6 Li, 7 Li, and 7 Be destruction, we found that 
the reaction channel which proceeds through the 1 + , 
E* = 0.770 ± 0.010 MeV nuclear excited state of 8 B via 
r Be x +P -> 8 B*(1+, 0.770 MeV)x ^ 8 B x + 7 could also 
destroy some 7 Bex, and that the destruction processes 
6 Lix (p, 3 He) 4 Hex, 7 Lix (p,a) 4 Rex might also be signifi- 
cant. 

Our calculated binding energies of the X" parti- 
cle in 7 Bex and 8 Bx are respectively 1.488 MeV and 
2.121 MeV. Adopting these values without any correc- 
tion to the energy levels of the nuclear excited states 
of 8 Bx, this 1 + state of 8 Bx is located near the parti- 
cle threshold for the 7 Bex+P separation channel. Thus, 
the 7 Bex(p,7) 8 Bx reaction can proceed through a zero- 
energy resonance of 8 B X . However, the measured energy 
uncertainty of the 1 + state of 8 B is ±10 keV, and more- 
over, the excitation energy of this level is very sensitive 
to the model parameters used to calculate the binding 
energies of the X-nuclei. Even such a small uncertainty 
of the resonance energy as 10-100 keV would dramat- 
ically change the TRR because the BBN catalyzed by 
the -nuclei proceeds at effective temperatures as low as 
Tg ~ 0.1. Taking account of the uncertainties associated 
with the 1 + resonance energy, E, from the 7 Bex~\-p sep- 
aration threshold, we found that E « 30 keV maximizes 
the TRR. This threshold energy would be achieved when, 
for example, the uniform charge radii are 2.2955 fm for 
7 Bex and 2.4564 fm for 8 Bx, respectively. This reso- 
nant reaction is potentially as effective as 7 Bex+P — > 
8 B^ ^ s Bx + 7 in destroying 7 Be. However, the charge 
radii we have adopted tend to be smaller than the mea- 
sured charge radii, and this might overestimate the bind- 
ing energies of X" . If a more realistic calculation were 
performed, the resulting binding energies might shift this 
1 + excited state upward, which would diminish the effect 
of this destruction process. In addition, the transition 
through this state would be E2 or Ml, which might also 
weaken its effect. Even in this case, though, the atomic 
resonance 8 B^f [8] plays the important role in destroying 



7 Be x . 

Since it is important to know precisely when during 
BBN the X~ particles become bound to nuclei, and what 
their distribution over the BBN nuclei would be Q at 
any time it is necessary to consider the thermodynamics 
associated with binding the X~ particles. We thus in- 
cluded both recombination and ionization processes for 
X~ particles in our BBN network code and dynamically 
solved the set of rate equations to find when the A-nuclei 
decoupled from the cosmic expansion. 

Regarding the thermonuclear reaction rates we note 
that since the mass of X~ particle mx is assumed to be 
> 50 GeV, the reduced mass for the X~+A(N, Z) system 
can be approximated as fix = mAmx / (mA+mx) ~ ttia, 
rendering the thermonuclear reaction rate for the first 
recombination process A(X~ ,j)Ax @| to be (a r v)x ~ 
2 9 7raZ 2 (2^) 1 / 2 /(3 exp(4.0))S bind /(^( M xT) 1 /2) 
m^ 2 ' 5 , where a is the fine structure constant. This 
rate is almost independent of mx- However, the rate 
for the second recombination process Ax{X~, j)Axx 
is dependent upon mx, i.e., (<7 r v)xx ~ 2 9 ira(Z — 
l) 2 (27r) 1 /2/(3exp(4.0))i? bind /( A1 2 , Jf (^r) 1/2 ) 
m x 2 ' 5 . This arises because fixx = mAxmx /{rriAX + 
mx) ~ mx/2. Since mx is assumed to be much larger 
than the mass of the light nuclei mx 3> niAi the rate 
for the second or higher-order recombination process is 
hindered. 



III. RESULTS 

The evolution of the BBN abundances when X~ par- 
ticles are included exhibits some particularly notable fea- 
tures. During the nucleosynthesis epoch, the abundances 
for 6 Li, 7 Li and 7 Be assume their normal BBN values un- 
til the temperature reaches Tg ~ 0.5 — 0.2. Below that 
temperature the X~ particles bind to the heaviest nu- 
clides, 7 Li and 7 Be. When the abundance ratio, Yx, of 
X~ particles to baryons is larger than 0.1 these nuclides 
are then partially destroyed by reactions that would have 
previously been inhibited by the Coulomb barrier. At 
around Tg = 0.1, the X~ particles are captured onto 4 He. 
Then a new round of X-nuclei nucleosynthesis occurs. In 
particular, the reaction 4 Bex(d,X~) 6 Li produces normal 
6 Li nuclei with an abundance which is orders of mag- 
nitude above that from standard BBN. An interesting 
feature is that the 6 Li formed in this way is not easily 
destroyed by the 6 Li(p,a) 3 He reaction, the dominant 6 Li 
destruction reaction in BBN, because the X~ transfer 
reaction restores the charge to 6 Li. Hence, the Coulomb 
barrier is too high at this temperature for its destruction 
resulting in a large 6 Li/ 7 Li abundance ratio. 

The final calculated abundances of the mass 6 and 7 
nuclides, however, depend strongly on the assumed X~ 
abundance. At high X" abundance levels, more than 
one X" capture can occur. Although the abundance of 
these multiple X~ bound particles is too small to sig- 
nificantly contribute to BBN, they nevertheless interact 
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readily since their Coulomb barriers are greatly reduced. 
This is especially true of charge-neutral A B\e X x- To clar- 
ify the nucleosynthesis yields, we have thus made a study 
in which the X~ abundance Yx was varied over a wide 
range. 

In Fig. 1 we show contours of an interesting region in 
the decay lifetime tx vs. Yx plane. Curves are drawn 
for constant lithium abundance relative to the observed 
value in MPHSs, i.e., d( 6 Li) = 6 Li Calc/6 Li Obs ^ qM 
curves) and d( 7 Li) = 7 Li Calc / 7 Li obs (dashed curves) for 
several values of the stellar depletion factor "d". The 
adopted abundances are 7 Li/H= (1.23±g;||) x 1(T 10 
and 6 Li/Fh= (7.1±0.7) x 1CT 12 1]. Shaded regions for the 
d( 6 Li) = 1 and d( 7 Li) = 1 curves illustrate the 1 a uncer- 
tainties in the adopted observational constraints based 
upon the dispersion of the observed plateaus. We also 
show curves for stellar depletion factors of d( 7 Li) = 2,3 
and d( 6 Li) = 4, 25. Since 6 Li is more fragile to stellar 
processing than 7 Li [l^|, its possible depletion factors 
could be larger than those for 7 Li. 

The main point of this figure is that, independent of 
stellar destruction, it is possible to find a simultaneous 
solution to both the 7 Li overproduction problem and the 
6 Li underproduction problem. This occurs in the param- 
eter region Y x « 0.09 - 0.6, t x « (1.6 - 2.8) x 10 3 s 
consistent with the suggestion of Q. Assuming that 
the products of the decaying X~ particles are progen- 
itors of the CDM particles, the WMAP-CMB obser- 
vational constraint on f^cDM = 0.2 limits the mass of 
the X~ , i.e., Y x m x < 4.5 GeV and m x < 50 GeV 
when we include the destruction reaction processes of 
A = 7 nuclide 7 Be x +P -> 8 Bf -^ 8 B X + 7 @ and (as- 
sumed maximal value of the) 7 Bex + p ^ 8 B*(1 + , 0.770 
MeV)x ^ 8 Bx+7- When we include the destruction pro- 
cess 7 Bejc ^ 7 Li+X° [8J] , these parameter ranges slightly 
change to Y x « 0.04 - 0.1, t x « (1.8 - 3.2) x 10 3 s, and 
m x ^5 100 GeV. Figure 2 illustrates the final calculated 
BBN yields as a function of baryon to photon ratio r\ 
for the case of (Yx, tx)=(0.6, 1.6 x 10 3 s). This choice 
leads to 6 Li and 7 Li abundances consistent with the ob- 
served values without stellar depletion. Note, though, 
that the same conclusion is reached if the destruction re- 
action through the 8 B(1 + ) state is not included, so the 
general conclusion is robust. 



IV. SUMMARY 

In summary, we have investigated light-element nucle- 
osynthesis during BBN taking into account the possibility 
of massive, negatively-charged X~ particles which would 
bind to the light-nuclei. When the chemical and kinetic 
processes associated with such particles are included in 
a BBN code in a fully dynamical manner, along with the 
reactions enabled by the X~ particles, the X~ particles 
are found to enhance the reaction rates in BBN, both 
by reducing the charge of the resulting A-nuclei, and 
by enabling transfer reactions of the X~ particles. X~ 
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FIG. 1: Contours of constant lithium abundance relative to 
the observed value in MPHSs, i.e., d( 6 Li) = 6Li Calc / 6 Li° bs 
(solid curves) and d( 7 Li) = 7 Li Calc / 7 Li° bs (dashed curves). 
The adopted abundances are 7 Li/H= (1.23±o™) x 10" 10 
and 6 Li/H= (7.1 ± 0.7) x 10" 12 [J. Shaded regions for the 
d( 6 Li) = 1 and d( 7 Li) = 1 curves illustrate the 1 a uncertain- 
ties in the adopted observational constraints based upon the 
dispersion of the observed plateaus. 

particles greatly enhance the production of 6 Li, primar- 
ily from the X~ transfer reaction 4 Hex(d,X~) 6 Li. The 
7 Li abundance, however, decreases when the X~ parti- 
cle abundance is larger than 0.1 times the total baryon 
abundance. In this case, the 7 Li abundance decreases 
with the X~ particle abundance due to the inclusion of 
two resonance channels for 7 Bej5f(p,7) 8 Bx through the 
nuclear and atomic excited states of 8 Bx- It was found 
to be important to predict precisely the binding energies 
and excited states of exotic A-nuclei in realistic quantum 
mechanical calculations. Both abundance ratios of 6 Li/H 
and 7 Li/H observed in MPHSs are obtained with an ap- 
propriate choice for the lifetime and abundance of the X~ 
particle. These observational constraints imply a lifetime 
and abundance roughly in the range of tx ~ 2 x 10 3 s and 
Yx ~ 0.1. We deduce that this Yx value requires that 
mx ~ 50 GeV in order to guarantee that this abundance 
of A~ particles survives to the epoch of nucleosynthesis. 
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FIG. 2: Abundances of 4 He (mass fraction), D, 3 He, 7 Li and 
6 Li (by number relative to H) as a function of the baryon-to- 
photon ratio rj or QbIi 2 . The dashed and solid curves are re- 
spectively the calculated results in the standard BBN and the 
X' catalyzed BBN for the case of (Y x , rx)=(0.6, 1.6 x 10 3 s). 
There is virtually no difference between the dashed and solid 
curves for 4 He, D, and 3 He. The band of theoretical curve for 
each nucleus displays 1 a limits taken from [l5l ]. The hatched 
regions represent the adopted abundance constraints from [Trj | 
for 4 He, for D, for 7 Li, and fl] for 6 Li, respectively. 
The vertical stripe represents the 1 a Qsh 2 limits provided 
by WMAP @]. 
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